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ABSTRACT 

Wc present a measurement of the angular power spectrum of the cosmic microwave background 
(CMB) radiation observed at 148 GHz. The measurement uses maps with 1.4' angular resolution 
made with data from the Atacama Cosmology Telescope (ACT). The observations cover 228 deg^ 
of the southern sky, in a 4? 2- wide strip centered on declination 53° South. The CMB at arcminute 
angular scales is particularly sensitive to the Silk damping scale, to the Sunyaev-Zel'dovich (SZ) effect 
from galaxy clusters, and to emission by radio sources and dusty galaxies. After masking the 108 
brightest point sources in our maps, we estimate the power spectrum between 600 < i < 8000 using 
the adaptive multi-taper method to minimize spectral leakage and maximize use of the full data set. 
Our absolute calibration is based on observations of Uranus. To verify the calibration and test the 
fidelity of our map at large angular scales, we cross-correlate the ACT map to the WMAP map and 
recover the WMAP power spectrum from 250 < I < 1150. The power beyond the Silk damping tail 
of the CMB ~ 5000) is consistent with models of the emission from point sources. We quantify the 
contribution of SZ clusters to the power spectrum by fitting to a model normalized to cts = 0-8. We 
constrain the model's amplitude Asz < 1-63 (95% CL). If interpreted as a measurement of erg, this 
implies erf ^ < 0.86 (95% CL) given our SZ model. A fit of ACT and WMAP five-year data jointly to 
a 6-parameter ACDM model plus point sources and the SZ effect is consistent with these results. 
Subject headings: cosmology: cosmic microwave background, cosmology: observations 
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1. INTRODUCTION 

The cosmic microwave background (CMB) radiation 
captures a view of the universe at only ~ 400, 000 years 
after the Big Bang. The angular power spectrum of 
temperature anisotropics in the CMB has been crucial 
in developing the current standard cosmological model, 
in which the universe today contains some 5% baryonic 
matter, 23% dark matter and 72% dark energy. We re- 
fer to this model throughout as the ACDM model. The 
temperature power spectrum has been measured to good 
precision for multipole moments £ < 3000. At angular 
scales with £ < 2000, the power spectrum matches the 
predictions of ACDM, and it can be used to constrain 
multiple parameter s of the cosmological model (e.g., 
Dunklev et all 120091: iB^own et al.l[200l iReichardt et alJ 



uunKiey et ai.M:^UUa Itirown et ai.M^UUa meicnarat et ai.l 
2009bt iSievers et al.l 120091 ). The agreement between 



the current polarization anisotropy measurements and 
the predictions of ACDM cos mology further supports 
the model (e g.. iDunklev et al.l [20091: iBrown et"aLll2009l: 
IChiang et al.ll2009D . 

At ^ > 3000, the signal from the primary CMB 
anisotropy becomes dominated mainly by two popu- 
lations. The first consists of point source emission 
from both radio and dusty infrared-emitting galaxies. 
The second is the population of massive galaxy clus- 
ters that give rise to the Sunya ev-Zel'dovich (SZ) ef- 
fect (jSunvaev fc Zel'dovichl ri970( ). in which CMB pho- 
tons scatter off the electrons of the hot intra-cluster 
medium. Removal of the brightest foreground galax- 
ies and SZ clusters can reduce their contributions to 
the power spectrum, but the net anisotropy power due 
to unidentified sources still dominates the exponentially 
falling primary anisotropy spectrum at small scales. 

The top panel of Figure [T] shows recent measurements 
of the anisotropy at £ > 2000. In analyses of 30 GHz 
data, radio point sources are masked out and a resid- 
ual component is modeled and subtracted. For 150 GHz 
analyses, radio sources are masked and the residual radio 
and dusty galaxy contribution is estimated. 

The SZ effect has a unique frequency signature. In 
CMB temperature units (the units of Figure [T]) the am- 
plitude at 150 GHz is roughly half that at 30 GHz. The 
amount of SZ power is governed by erg, which measures 
the amplitude of the cosmic matter power spectrum on 
8/).~^Mpc scales; the SZ power scales approximately a s 
alinth)^ (ISeliak et al.l [2001 iKomatsu fc SeliakI 120021) . 
At 30 GHz, CBI reports ([Si^rs et al.ll2009n ' excess emis- 
sion above the ACDM model at £ » 3000 after acco unt- 
ing for all known radio sources ason "eFaI[[2009l) . A 
possible source of the "CBI exce ss" is the SZ effect. 
The SZA data ([Sharp et al.l [2001 ). also at 30 GHz and 
also after accounting for point sources, are consistent 
with the ACDM model at £ = 4000. The ACBAR re- 
sults ([Reichardt et al.ll2009bl) at 15 GHz are consistent 
with the CBI excess and ACDM. ISievers et all (|2009[ ) 
show that all the above data in the top panel of Figure [T] 
are consistent, within 95% CL, with ACDM plus a SZ 
contribution of erg = 0.922 ± 0.047 {la error bars). Re- 
cently, the South Pole Telescope (SPT) group reported 
(jg = 0.773 ± 0.025 based on the power spectrum at 150 
and 220 GHz ([Lueker et al.|[2009t) . 

In this paper, we present a new measurement of the 
CMB anisotropy power spectrum in the range 600 < £ < 



8000 (corresponding to angular scales of approximately 
1.4' to 18'). The observations were made at 148 GHz in 
2008 with the Atacama Cosmology Telescope. Figure [T] 
(bottom panel) shows the results, which are discussed in 
Section [7| We briefiy describe the instrument, the obser- 
vations, the calibration of the data, and the map-making 
procedure. The dynamic range of the power spectrum 
is large enough that we employ new techniques in spec- 
tral estimation. We confirm that difference maps of the 
data give power spectra consistent with no signal and 
that the power spectrum is insensitive to several details 
of our analysis. We use the shape of the power spectrum 
to bound the dusty galaxy contribution and the power 
from SZ clusters. 

2. INSTRUMENT AND OBSERVATIONS 

The Atacama Cosmology Telescope (ACT) is 
a 6-meter, off-axis Gregorian telescope optimized 
for arcminute-scale CMB anisotrop y measurements 



([Fowler et al.ll2007t iHincks et afl 120091 ). It was installed 
at an elevation of 5190 m on Cerro ToccQ in the Atacama 
Desert of northern Chile in March 2007. Observing 
conditions in the Atacama are excellent owing to the 
elevation, the arid climate, and the stability of the at- 
mosphere. After all cuts, the median precipitable water 
vapor (PWV) was 0.49 mm during the observations 
presented here. 

The MiUimeter Bolometer Array Camera (MB AC), 
the current focal-plane instrument for ACT, uses high- 
purity silicon lenses to reimage sections of the Grego- 
rian focal plane onto three rectangular arrays of detec- 
tors. The arrays each contain 1000 transition edge sensor 
(TES) bolometers. Th eir spectral cove rage is determined 
by metal-mesh filters ([Ade et al.ll2006f ) having measured 
band centers of 148 GHz, 218 GHz, and 277 GHz. The 
bolometers are cooled to 300 mK by a two-stage helium 
sorption fridge backed by commercial pulse-tube cry- 
ocool ers. The telescope performance and control sys- 
tems |Hincks et al."2008VSwitzer et al."2008'), camera de- 
sign (^hornton et al. 2008; Swetz et al. 2008), detector 
properties (iZhao et al. | |2008; Niemack et al. 2008), and 
readout electronics ([Battistelli et al.,,2008) are described 
elsewhereO 

2.1. Observations 

ACT operated with all three arrays from mid- August 
to late December, 2008 for the data presented here. The 
observing time was divided between two regions away 
from the galactic plane. The deepest observations cover 
900 deg^ of the southern sky in a strip 8° wide centered 
on (5 = -53° with RA from 19*^ to 24'' and O'^ to 7''36"\ 
The analysis presented in this paper uses only data from 
the central 228 deg^ of the southern strip and only ob- 
servations made with the 148 GHz array. The slightly el- 
liptical beam has full- widths at hal f-maximum (FWHM ) 
of 1.40' by 1.34' at this frequency ([Hincks et al.ll2OO90 . 

The observations were made by scanning the sky at a 
constant elevation of 50°. Each scan is 4? 5 wide on the 
sky (7?0 in azimuth angle). Scans repeat every 10.2s. 
Each half-scan consists of 4.2 s of motion at a constant 

1 ACT is at 22.9586° south latitude, 67.7875° west lon gitude. 
^ The site Ihttp : //www. physic s .prlnceton. edu/ act / 1 archives 
papers by the AC'i' collaboration. 
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Fig . 1. — Recent measurements of the CMS power spectrum, including this w ork. Top: the measurements of WMAP fNolta et ahl 
[200gt). Bolocam (Savers ct al. 2009), QUaD (Brown et al. 2009; Friedman et al. 20091). APE X-SZ (Rcicliardt et al. 2009a), ACBAR 
HReichardt et al.ll2009h) . SZA dSharp ct al.ii200al . BIMA (iPawson et al..2006i ). CBI HSievers et al.t200a) . and SPT (iLueker et al.,.200a ). For 
all the results, a radio point source contribution has been removed either by masking before computing the power spectrum (at 150 GHz), 
or by masking and modeling the residual (at 30 GHz and for WMAP). APEX-SZ additionally masks clusters and potential IR sources. 
Bottom: The ACT power spec trum from this work. The inset shows the cross-power spectrum between ACT and WMAP maps in the ACT 
southern field (sec Section 13.51 1 , which we use to check both the validity of the maps at larger scales and the absolute calibr ation. Only the 
ACT power spectrum is analyzed in this paper. In both panels and the inset, the solid curve (blue) is the ACDM model of lDunklcv et ahl 
II2009I ) (including lensing). The SZ effect and foreground sources are expected to contribute additional power, as shown in Figure |4] and 
Table [T] For display purposes — and only in this figure — we scale our result by 0.96 in temperature relative to the Uranus calibration; this 
calibration factor best fits our data to the ACDM model and diffe rs fr om the Uranus calibration by O.Tcr. Recent WMAP observations of 
Uranus suggest the same rescaling factor (see footnote to Section 13. 2( . ACT bandpowers for £ > 4200 have been combined into bins of 
A£ = 600 for this figure; they are given in a note to Table [T] 



speed of 1?5 /s followed by 0.9 s of acceleration. The first 
half of each night is spent observing the field rising in the 
eastern sky, after which ACT turns to the western sky 
to observe the same field as it sets through the standard 
elevation of 50°. The scan strategy is designed to mini- 
mize changes in the telescope's orientation with respect 
to the local environment while ensuring cross-linked ob- 
servations in celestial coordinates. Sky rotation ensures 
that all detectors sample all points in the field each night, 



apart from small areas at the edges. 

As the telescope scans in azimuth at constant elevation, 
each detector is sampled at 399 Hz. The data sampling, 
position reading, and all housekeeping data arc synchro- 
nized by a shared 50 MHz clock; absolute times arc ref- 
erenced to a GPS receiver with 0.25 ms accuracy. The 
data arc stored in continuous fiftccn-minutc segments 
called time-ordered data sets (TODs). Each TOD re- 
quires 1.6 GB of storage per detector array, or 600 MB 
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after applying lossless compression. 

In addition to the main survey, we perform occasional 
calibration measurements. Most nights, a few minutes 
are used to measure a planet when it passes through our 
standard elevation. Any given target is used only once 
in every three nights to minimize non-uniformity in the 
coverage of the CMB regions. The planetary observa- 
tions allow us to measure the system's relative and abso- 
lute responsivities, pointing, and beam profiles; they also 
provide a way to check for any time variations in each. A 
series of tuning and biasing procedures are also followed 
each night before and after regular CMB observations, 
to optimize and roughly calibrate the detector response. 

3. DATA REDUCTION AND MAP-MAKING 

The goal of the data reduction is to estimate the 
maximum-likelihood map of the sky. We select properly 
tuned detectors and calibrate their pointing and relative 
gains. The absolute gains are determined by observa- 
tions of Uranus. The map-making itself is an iterative 
and computationally intensive process. 

3.1. Data selection 

The first step in the data reduction is to select TODs 
when the receiver is cold enough for stable operation and 
when the precipitablc water vapor (PWV) is less than 
3 mm. Of the « 2880 h between Aug 25 and Dec 24, 
2008, MBAC was on line for 1352 h, roughly correspond- 
ing to the night time hours. After cutting on PWV and 
instrument performance, 1031 h (35% of calendar time) 
went into the pipeline. Of those, 850 h were spent ob- 
serving the southern strip. 

In the second step of the data reduction, we check 
whether each detector's data exhibit problems that war- 
rant removal from the final analysis (|Dunne 3|200l). In 
the case of rare and transient effects, data are removed 
for a single detector up to a few seconds at a time; other 
problems are diagnosed in 15-minutc intervals, and when 
necessary, each detector's data may be removed for the 
entire TOD. We check the detector's response speed; look 
for transient effects such as cosmic ray hits; and en- 
sure that the feedback loop remained loc ked, keeping the 
SQU ID amplifiers in the linear regime (jBattistelli et al.l 
|2008[) . Finally, we compare the detector's response to 
atmospheric emission with that of the array-wide aver- 
age, in each case filtering out data at frequencies above 
50mIIz. Detectors are cut if their data do not have at 
least a 0.98 correlation with the average or if the ampli- 
tude of the atmospheric signal is not within 15% of the 
median amplitude found among good detectors. This 
last cut is the largest and typically removes 100 detec- 
tors. Additionally, we bin the data by scan and compute 
the rms, skewness, and kurtosis in each bin. A cut on 
rms removes an additional 30 detectors on average, while 
the combined cut on skewness and kurtosis removes only 
one additional detector. When the instrument is well- 
tuned, ACT typically has 680 science-grade detectors at 
150 GHz. The remainder are discarded from the analy- 
sis or used to monitor instrumental effects. The effective 
sensitivity of the array, including all instrumental and at- 
mospheric effects, is ^ 30 |J.K s^/^ in CMB temperature 
units. 

3.2. Calibration to planets 



After we have a reliable set of detectors to use, the 
next step is gain calibration. In practice, the cuts are 
computed twice, so that the final set of cuts is based on 
calibrated data. We use Uranus as our absolute reference 
standard. 

The atmosphere is an excellent continuous fiat-fielding 
calibrator. We find that the time variation in detector 
relative responsivities is not significant. Therefore, the 
detectors' relative responses to the atmosphere (averaged 
over the entire season) are used to convert each output to 
a common scale. Observations of Saturn confirm these 
relative calibrations; it is the only source observed in 
2008 bright enough to use for this confirmation. 

We take the brightness temp erature of Uranus at 
148 GHz to be = 112 ± 6 K (iGriffin fc OrtonI [l99l 
iMarten et 311120051 iKramer et al.ll2008[ ). All results in 
this work scale in proportion to this temperature. For 
planet calib rations, we use the beam solid angle ^Ia — 
218 ± 4nsr (jHincks et al.ll2009[) . The net calibration er- 
ror of 6% in temperature is the combination of these 
two uncertainties and is dominated by uncertainty in the 
brightness of Uranus. 

Saturn was used as a rough check on the absolute cal- 
ibration result. The Saturn observations suggest that 
the ACT data are 10% brighter than inferred from the 
Uranus data. We attribute this inconsistency to the dif- 
ficulty in modeling the brightness of the planet and its 
rings as they vary over time and across frequencies. 

We track the stability of the system over time in a num- 
ber of ways. The absolute celestial calibration is checked 
through measurements of a planet on most nights. The 
conversion between raw data and units of power absorbed 
on the detectors is calibrated in two ways. The conver- 
sion is estimated twice nightly by modulating the TES 
detector bias voltage with a small additional square- wave 
and measuring the response. A calibration is also per- 
formed at the beginning of each night by sweeping all 
detectors through the full range o f bias vol t ages. Both 
methods are described further in [Niemack' (2008|) and 
IHshcr (200i). Based on several tests (Switzcr 200i), we 
find that the gain of each detector is constant over the 
season to better than 2%. 

As part of the calibration, the detector temporal re- 
sponse is deconvolved from the detector time streams. 
We model the response with a single time constant. The 
detectors' time constants are determined from the planet 
measurements and bias steps (Hincks 2009) . The median 
time constant is Tdet = 1.9ms (thus /sdB = 84 Hz), which 
corresponds at the ACT scan rate to £ ^ 31,000. Slow 
detectors with /sdB < 15 Hz are not used. We also decon- 
volve the anti-aliasing filter imposed in the data acqui- 
sition system ([Battistelli et al. 2008). We then further 
filter and sample the time stream at half the raw rate to 
speed the map-making step. This last anti-aliasing filter 
acts only at angular scales smaller than the beam and 
does not affect the maps. 

3.3. Pointing reconstruction 

The planet Saturn is bright enough that each bolome- 
ter detects it with high signal-to-noise in a single 
scan. We fit in the time domain to find the best two- 
dimensional location for each detector relative to the 
notional array center. The fit produces altitude and 
azimuth offsets for each, along with information about 
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Fig. 2. — The map and a difference map of the ACT southern field at 150 GHz. The same filters used in the power spectrum analysis 
are applied to both maps: an isotropic high-pass filter suppresses power for £ < 300 (Equation [2J , and all modes with \£x\ < 270 are set to 
zero, as described in Section 15.21 Top: The ACT southern field. The intensity scale is yiK (CMB units). The 360 deg^ with lowest noise 
are shown. The squares 4? 2 on a side indicate the 13 patches used for the CMB power spectrum analysis (228 deg^ total). Middle: A 
difference map made from two halves of the same data set. Most of the remaining structure visible at large scales is well below the range 
of £ that we consider in the power spectrum analysis. Bottom: The rms temperature uncertainty for one square arcminute pixels. 



the relative gains, the beam sizes, and the detector 
time constants. The detector pointings are consistent 
with optical models of the telescope and reimaging op- 
tics (iFowler et al.|[2007| ). The relative pointing used in 
this work is an average over 22 observations of Saturn. 
The rms uncertainty is 1.2" for the detector relative 
pointings, which contributes negligibly to any pointing 
error in the sky maps. 

The location of the array center can also be found from 
each observation of Saturn. It has a scatter of 4.3" rms 
over the 22 observations, which we attribute to slight 
thermal deformations in the telescope mirrors and their 
support structures. 

The Saturn data give the absolute location of the ar- 
ray center only for other observations taken at the same 
horizon coordinates. To determine the absolute location 
of the array center during science observations, we used 
approximately 20 known radio sources found in prelimi- 
nary maps. We find this approach both simpler and bet- 
ter constrained than making a complete pointing model 
of the telescope. The ACT southern field was observed 
both rising and setting, at azimuths centered at 30° on 
either side of south. These preliminary maps were made 
separately for rising and setting data and were used to 
determine separate pointing corrections for the two cases. 
After correction, we estimate the maps to have 5" point- 
ing uncertainty. 

3.4. Map-making 



The goal of the map-making step is to take the 3200 GB 
of cut, calibrated, and deconvolved raw data from both 
the rising and setting scans and produce from it a maxi- 
mum likelihood estimate of the sky. This work produces 
a 200 MB map, 16,000 times smaller than the raw data 
set. We first multiply the data for each detector in each 
TOD by a window function, reducing the weight in the 
first and last 10 seconds of each file, then remove a single 
offset and slope for the entire 15-minute period. There 
is no additional filtering in the time domain, though the 
lowest frequencies are given no weight in the process of 
maximizing the likelihood of the final map. The detector 
data are combined into a data vector d. The mapping is 
done in a cylindrical equal-area projection with a stan- 
dard latitude of J = -53.5° and pixels of 30" x 30", 
roughly one third of the beam FWHM. The map is rep- 
resented as a vector m of length A'pix ~ 10^. We model 
the data as d = Pm 4- n where the matrix P projects 
the map into the time stream, and n is the noise, which 
has covariance matrix N. The maximum likelihood solu- 
tion, m, is giv en by solving the mapping equation (e.g., 
lTegmarkl[T997i) : 

P^N-iPm = P'^N^M. (1) 

We solve for rh iteratively using a preconditioned 
conjugate gradient (PCG) method ([Press et al.l 120071 : 
iHinshaw et al.l I2007D . Based on simulations, we find 
that the solution is an unbiased estimator of the sky for 
£ > 600, the multipoles that we analyze in this paper. 
Additionally, the clear cross-correlation with WMAP 
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(Section I3.5|) indicates that the maps are likely unbi- 
ased as low as £ ^ 200. This approach is different 
from that taken by other recent measurements of the 
fine angular scale anisotropy in which the time stream 
data are fi ltered and binne d fe.g.. iFriedrnan et al.ll2009l : 
IReichardt et al. 2009a: Lueker et al.ll2009[) . 

Solving Equation[T]requires careful consideration of the 
noise structure of the ACT data. First, the data are 
weighted in Fourier space to account for the variation in 
noise with frequency, particularly to reduce the noise that 
the atmosphere adds at low frequencies. For frequencies 
less than 0.25 Hz, we set the statistical weights to zero. 
Higher frequencies are given successively more weight in 
inverse proportion to the noise variance in each band. 
Second, we find that there are several modes, particular 
combinations of the 1000 detectors, that correspond to 
signals unrelated to the celestial temperature; for exam- 
ple the common mode, atmospheric gradients, and other 
detector correlations. We handle them by finding the 
ten modes with the largest eigenvalues in each TOD and 
solving for their amplitudes as a function of time along 
with the map rh. The power spectrum (Section [5]) is 
robust to halving or doubling the number of modes. 

Maps are made using only fractions of the valid data 
so that null tests can be performed. The subset maps 
are also used for finding th e CM B power spectrum, as 
described further in Section [5^ The final map used for 
source-finding, however, is a complete run based on all 
the valid data. Figure [5] shows the map and the region 
used. The computational task is considerable. One iter- 
ation takes 100 seconds on 5000 cores of Canada's SciNet 
GPC cluster each running at 2.53 GHz; computations re- 
quired before the first iteration take an additional time 
equivalent to approximately twenty iterations. A con- 
verged map requires hundreds of iterations and approxi- 
mately ten CPU years. 

3.5. Calibration to WMAP 

The absolute calibration described in Section 13.21 is 
limited by the uncertainty in the 148 GHz brightness 
of Uranus, the primary calibrator. We complement the 
planetary calibration by cross-correlatin g with the W- 
band (94 GHz) measurements of WMAP (jHinshaw et al.l 
120091) . In the angular scale range of 200 < ^ < 1200, 
both instruments measure the sky with sufficient signal- 
to-noise to permit the comparison. 

The technique requires a single map from WMAP and 
two maps with independent noise from ACT. We assume 
that only one relative calibration ratio between ACT and 
WMAP must be estimated. The WMAP data set con- 
sists of the high-resolutior0 W-band five-year map. The 
two independent ACT maps are each made with one half 
the current data set. We compare the maps only in the 
228 deg^ strip used for th e pr esent power spectrum anal- 
ysis (described in Section [57l|) . where the noise is low and 
uniform. 

We compute those cross-spectra in two-dimensional i- 
space that combine either the full ACT data set with 
the WMAP map or the two ACT maps with each other. 
We average over the polar angle in ^-space to get a one- 
dimensional cross-spectrum for each. In the case of the 



ACT spectra, the noise is not isotropic, so the spectra 
are angle-averaged with appropriate weights. We then 
find the single calibration factor that minimizes for 
a model in which the minimum-variance weighted com- 
bination of the ACT- ACT and the ACT- WMAP cross- 
spectra equals the WMAP all-sky spectrum. The model 
accounts for the fact that the ACT and WMAP mea- 
surements have noise varying as different functions of L 
The result is a calibration factor with < 6% fractional 
uncertainty in temperature. The calibration factors de- 
rived from WMAP and from Uranus observations are 
consistent to < 6%. 

4. FOREGROUNDS 

Radio and infrared galaxies are the dominant sources 
of foreground emission in the 150 GHz band &t t > 
1000. To study the underlying CMB spectrum, the 
sources must be identified and masked, and resid- 
ual contamination in the power spectrum must be 
accounted for. Many approaches to the problem 
have been described in the CMB lit er ature (e.g., 
Wright et al."2009; Rcichard t et~ani20 09a'.'b': 'S harp et alJ 
1)09; Dawson et al. 2006; Sievers et a l. 2009). Our ap- 
proach is to find the sources in the ACT maps, mask 
them, and assess the residual contribution with models. 

We find the sources using a matched filter with noise 
weightin g derived from the s tatistical properties of the 
map (Tegmark fc de Oliveira -Costa 1998). We identify 
as sources all pixels detected at > ba having at least 
3 neighboring pixels detected at > 3cr. A 5cr detection 
corresponds to roughly 20 mJy in the filtered map. The 
selection criteria are tested on simulated sky maps to 
determine the sample's purity and completeness. The 
sample is approximately 85% complete at 20 mJy and 
~ 100% complete at 50mJy. For fluxes greater than 
20 mJy, the detections have a purity of ^ 95% (that is, in 
a hundred detections approximately five are false) . Sim- 
ulations show that at 20 mJy deboosting (lCondonlll974D 
is a '--^ 9% effect, and at 40 mJy a < 1% effect. 

In the 228 deg^ area, we detect 108 sources. Of 
these, 105 can be ide ntified wi th sources in th e 
PMN (iWright et al.l[l99l. SUMSS (iMauch et al.l [2003h. 



3 The WMAP maps are at HEALPix resolution A^sido 
with 3.5' pixels. 



1024, 



and /or AT20G (jMurphv et al.|[20"09[) radio source cata- 
logs. One source, not in these catalogs, is in the 2MASS 
catalog (and thus has IR emission). Two have no previ- 
ously measured counterparts. In a 1' search radius, 17 
sources have both radio and IR identifications. With the 
exception of the single 2MASS source, all are treated as 
radio sources. 

We fit the fiux distribution of the radio sources to three 
model s. The distribution agrees with the lToffolatti et ahl 
(|1998| ) radio model scaled by 0.49 ±0.12 in source counts, 
where the uncertainty is statistic al only. In the WMAP 
analysis of sources at 41 GHz, iHinshaw etaD ()2007l ) 
found a good fit with the same model after scaling the 
counts by 0.64. Thus it appears that on average the ra- 
dio population at 150 GHz is composed predominantly 
of fiat-spectrum sources. We prefer to compare the data 
to models by scaling the model's flux rather than the 
counts, as the extrapolation of the radio fiux into the 
150 GHz range is ty pically the most uncer tain part of 
any model. For the iToffolatti et al.l (I1998D model, the 
requir ed flux scaling i s 0.5. The sources are also well fit 
by the lDe Zotti et all (|2005j ) model scaled by 1.2 in fiux. 
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Lastly , we fit the radio sources model in iSehgal et alj 
()2010j) and find that its flux should be scaled by approx- 
imately 2.1. All scale factors have a statistical uncer- 
tainty of about 25%. Each of these models can be used 
to estimate the residual rms flux below the 20mJy cut, 
as discussed in Section 17.11 We plan a more thorough 
exploration of the radio sources in an upcoming publica- 
tion. 

Before computing the power spectrum, we mask a 10' 
diameter region around each of the 108 sources. We call 
this Mask-ACT. In total, 2.3 deg^ are masked, or 1.0% 
of the map. Doubling the size of the masked holes has a 
negligible effect on the power spectrum. We also create 
a cluster mask, Mask-C, by finding all clusters and 
check that using this mask combined with Mask-ACT 
has a negligible effect on the power spectrum. We also 
generate a mask, Mask-R, based on known radio sources 
measured in the SUMSS (0.8 GHz), PMN (5 GHz) and 
ATCA (20 GHz) c atalogs. A niodel t hat has been fit to 
multiple data sets (jSehgal et al.ll20ldl ) is used to estimate 
the mean spectral index, and this index is then used to 
determine the flux cut in the radio catalogs. Since the 
model is based on ensemble properties of radio sources it 
may miss lower-flux low frequency sources having shallow 
indices. The resulting power spectrum is independent of 
which radio source mask we use, or if we use the union 
of the two, and does not depend on whether we include 
the cluster mask. 

Diffuse dust emission becomes dominant near the 
galactic plane. We check for evidence of dust in the 
region of our map nearest to the plane, from right as- 
cension 6^ to 7^. Through a cross-correlation with 
the estimated dust map of iFinkbeiner et all (fT999t ) (the 
FDS map), we can limit the dust contribution here to 
i{e + l)Q/(27r) < 5 at £ = 1000, while a direct 
power spectrum of the FDS map in the same region gives 
~ 1 |i.K^. As the diffuse dust component decreases with 
increasing £, it is not significant in our analysis. 

5. POWER SPECTRUM METHOD 

We estimate the CMB power spectrum using the adap- 
tive multi-ta per method (AM TM) with prewhitening, 
described in iDas et ahl ()2009( ) . We make independent 
maps from subsets of the data and use only cross-spectra 
between maps to estimate the final power spectrum. All 
operations are performed using the flat-sky approxima- 
tion. We summarize the method in this section. 

5.1. Fields used for power spectrum analysis 

We find the power spectrum of our map by separate 
analysis of each of the 13 patches shown in Figure [21 
Each patch is 4? 2 x 4? 2 in size, and together they cover 
a rectangular area of the map from a ~ 0''48'" to 6^52™ 
(12° to 103°) in right ascension and from (5 = —55° 16' 
to — 51°05' in declination. This area is the region of a 
larger survey having the lowest noise. We also split the 
raw data into four subsets of roughly equal size, with 
the data distributed so that any four successive nights 
go into different subsets. The four independent maps 
generated from these subsets cover the same area and 
have approximately the same depth. All maps are fully 
cross-linked. That is, they all contain data taken with 
the sky both rising and setting. 



5.2. Spectrum of a single patch 

We estimate the spectra of the 13 patches indepen- 
dently, before taking a weighted average to find the final 
spectrum. There are four independent sky maps, from 
which six cross-spectra and four auto-spectra are evalu- 
ated on each patch. We use a weighted mean of only the 
cross-spectra for the final spectral estimate. The weights 
depend on both the cross- and auto-spectra, as discussed 
below. 

Before separating the four maps into 13 patches, each 
beam-convolved map, Tf,(0), is initially filtered in Fourier 
space with a high-pass function Fc{(.). This filter sup- 
presses modes at large scales that are largely uncon- 
strained. These modes arise from a combination of in- 
strument properties, scan strategy, and atmospheric con- 
tamination. We choose a squared sine filter, given in 
Fourier space by the smooth function 





sin' 
1 



£ <C ^min 

^min *\ ^ ^max 



(2) 



where x{£) = (7r/2)(£ - 4iin)/(4iax - ^min)- We choose 
^min = 100 and £max — 500. The 13 patches of the map 
are treated separately from this point. The map of each 
patch is prewhitened using a local, real-space operation 
to reduce the dy namic range of its Fourier components 
(jDas et al.ll2009() . The prewhitening operation involves 
adding a fraction (2%) of the map to an approximation 
of its Laplacian. The Laplacian is computed in real space 
by taking the difference between the map convolved with 
disks of radius 1' and 3'. Then the maps are multiplied 
by the point source mask. The prewhitening step greatly 
reduces the leakage of power from low to high multipoles 
caused by the action of the point source mask on the 
highly colored CMB power spectrum. 

For each patch, we compute the six 2D cross-spectra 
and four 2D auto-spectra. The axes correspond to right 
ascension and declination. Each spectrum is computed 
using the adaptively weighted multi-taper method, using 
A^ta p = 5^ tapers ha ving resolution parameter A^ros = 3 
fsee lDas et al.|[2009() . Windowing the maps with 25 or- 
thogonal taper functions allows us to extract most of the 
statistical power available in the maps, at the expense of 
broadening the resolution in angular frequency by a fac- 
tor of approximately N^cs- At this stage, each 2D cross- 
spectrum, C^"^ , between submaps a and /3 on patch i 
incorporates the effects of the filter, prewhitening, ta- 
pering, the point source mask, and the beam; they are 
analogous to a "pseudo powe r spectrum" of an apodized 
map (e.g.. lHivon et al1l2002D . Each 2D spectrum is then 
averaged in annuli with a narrow range of |£| to give the 
binned pseudo-spectrum Bb, with 



10.(3 



PbfB, 



■iaf3 



(3) 



Throughout this paper we define Be = £{£+l)Ce/2TT. The 
binning function Pf,t, is set to one for pixels lying in an 
annular bin (indexed by b) of width A£ ~ 300 centered on 
1^1 = £b, and zero elsewhere. The size of the patches and 
the resolution of the tapers dictates the width of the bins. 
For our square patches of side s = 4?2, the fundamental 
frequency resolution in Fourier space is (5£ = 27r/s w 90. 
The application of the tapers with a resolution param- 
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eter TVVes = 3 further degrades the resolution to Nj-^sSi, 
so bins chosen to be smaller than A£ ^ 270 would be 
unavoidably correlated. The binning function is also set 
to zero where \£x\ < 270. This region of Fourier space is 
particularly sensitive to scan-synchronous effects, either 
fixed to the ground or in phase with the azimuth scan. 

The binned pscudo spectrum is then deconvolved 
with the mode-mode coupling matrix, which takes into 
account the combined effects of tapering and masking 
and can be computed exactly. Lastly, we divide by the 
£-space representations of the prewhitening filter, of the 
high- pass filter Fc{£b), and of the beam to obtain an un- 
biased estimate of the true underlying spectrum . 

This procedure has been tested with simulations. The 
number and resolution of the tapers are chosen as the op- 
timal balance between maximizing information and min- 
imizing bias caused by leakage of power. The simulations 
confirm that increasing the number of tapers beyond 5^ 
has a negligible effect on the spectrum errors. 

5.3. Combining patches 

The final power spectrum estimator is given by a 
weighted mean over iV = f 3 patches. 



(4) 



wher^l Bl = J2a /3-a<p ^l"'^ mean of the six 
deconvolved cross-power spectra in patch i (assuming 
equal weights), and a and /3 index the four independent 
maps of that patch. The weights are chosen as the in- 
verse of the variance of this estimator in each patch, i.e. 
wl = l/a'^{Bl), where 



cr'iBl 



i\2 



Bl 



(5) 



and the average is taken over the several cross-spectra 
computed for patch i. The first term contains the 4-point 
function of the temperature field and is approximated as. 



M{M - 1) 



X {BfBl' + BfB^,^ + BrB^,'),{Q) 

where M = 4 is the number of submaps per patch. 
We have neglected any connected (non-Gaussian) part 
of the 4-point function due to components such as point 
sources. This is a reasonable approximation when choos- 
ing the 13 weights, because the expression (Equation [5]) 
is dominated by the auto-spectrum terms, which in turn 
are noise-dominated. 

5.4. Power spectrum covariance 
We estimate the bandpower covariance matrix S using 



the scatter in the power spectrum among the 
patches. 



^bb' 



- {ABbABb 



13 



(7) 



* Here we introduce the notation X as an unbiased estimator of 
the quantity X, in the sense (X) = X. 
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— Analytic total 

— Measurement noise (N) 

— Sample error (S) 

— Poisson error (P) 
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Fig. 3. — The total estimated error a on the power spe ctrum 
(blue) given by the analytic expression (Equations [9] and (TO]!. The 
uncertainty on cr is found using simulations and is shown by the 
shaded blue band. The Gaussian sample variance (red line, labeled 
"S") dominates for £ < 1200, and atmospheric plus instrument 
noise (purple, "N") dominates at £ > 2000. The non-Gaussian term 
due to unmasked point sources and clusters of galaxies contributes 
about 15% of the variance at 2500 < £ < 6000 (green, "P"). The 
errors estimated using the scatter of the results from the thirteen 
patches (Equation |8] are shown for comparison (black points); they 
agree well with the analytic errors. 



b'^b' 



{Bl-Bb){Bl,-Bb') 



N 



N 



(8) 



wl 



The square roots of the diagonal elements of the covari- 
ance matrix are reported as the errors on our power spec- 
trum estimate. 

To test the accuracy of the error estimate, the errors 
on the power spectrum are also computed analytically. 
Three terms contribute: sample variance in the CMB 
multipoles due to limited sky coverage, instrumental and 
atmospheric noise, and a non-Gaussian term due to un- 
masked point sources and galaxy clusters. The diagonal 
component of the variance in one patch can be written 
as the sum of these terms, in order: 



a [Lb) = 

Ub 



Ub 



LP 

/sky 



^, (9) 



where Nb — {C^ — Cb), the difference between the 
weighted mean auto- and cross-spectra, estimates the av- 
erage power spectrum of the noise; M = 4 is the num- 
ber of maps with independent noise properties per patch; 

= M{M— 1)/2 = 6 is the number of cross-spectra per 
patch; Ub counts the number of Fourier modes measured 
in bin h (that is, the number of pixels falling in the appro- 
priate annulus of Fourier space); /sky is the patch area 
divided by the full-sky solid angle. An steradians; and Cb 
is the weighted mean cross-spectrum. In the last term, 
(Tp is given by the non-Gaussian part of the four-point 
function. Thus 



1 

47r 



{5T\h))-2,{6T\h)y 



(10) 



where T is the temperature map and the average is over 



10 



J. Fowler et al. 



pixels in the map. The term arises from the Poisson- 
distributed components in the maps, including unre- 
solved point sources and clusters of galaxies. For pur- 
poses of the covariancc calculation, we assume that such 
sources are distributed independently of one other. This 
term is constant as a function of £ and is computed from 
the masked maps with the high-pass filter (Equation [2]) 
applied. The overall variance is given by the weighted 
mean of the patch variances (Equation |9l). The covari- 
ancc among bins is small in the limit that the four-point 
term does not dominate, and with bins chosen to be suf- 
ficiently large. 

Figure[3]shows the three components of the errors. The 
Gaussian part due to sample variance dominates below 
£ = 2000; it scales as Cb/y^. The Poisson term due to 
the non-Gaussian clusters and unresolved point sources 
contributes about 15% of the error budget for 2500 < 
£ < 6000 but is sub-dominant at all scales. The noise 
term dominates at ^ > 2000; atmospheric noise is the 
main contribution at angular scales of ^ < 1000. As a 
cross-check of the error estimates, we also plot the errors 
derived from the scatter of the power spectra from the 
thirteen patches. 

The uncertainty in the analytic errors is shown as the 
shaded band in Figure [3] It was estimated by Monte 
Garlo simulations. One thousand patches were simu- 
lated with white noise and their power spectra taken by 
the same methods used on the ACT maps. The results 
demonstrate that the errors estimated from the scatter 
among our 13 patches are consistent with the expected 
uncertainty. The same simulations were used to verify 
that the covariance between different power spectrum 
bins (Equation[8]) is less than 1%. 

6. POWER SPECTRUM RESULTS 

The binned estimate of the power spectrum Bb is shown 
in Figures 2] and [1] and bandpowers are given in Table [TJ 
With our method the bandpowers are estimated to have 
less than 1% correlation with neighboring bins, but the 
window functions have a small overlap, which we account 
for in the analysis. At ^ < 2500, the estimated power 
is consistent with previous observations by ground and 
balloon-based experiments. The features of the acoustic 
peaks are not distinguished with the coarse binning, but 
with the fluctuation band-power measured to 5%, this 
spectrum offers a powerful probe of cosmological fluctu- 
ations at small scales. A clear excess of power is seen at 
£ > 2500 which can be attributed to point sources and, 
to a lesser extent, to the SZ effect. 

The large-scale modes are recovered only after iterat- 
ing the maps. To test for convergence of noiseless maps, 
we compute the spectra CI of difference maps, between 
the processed simulation map and input simulated sky, at 
successive iteration numbers i. To account for differences 
that are non-uniform between data subsets, we use the 
auto-spectra. We show that the amplitude of fluctuations 

in the difference map {^J^l) is small, less than 1% of the 

amplitude of fluctuations in the input map (VC^) at all 
scales by iteration 500, as shown in Figure [SJ To test 
for convergence in the data, where we do not know the 
input map, we estimate the maximum change in power 
between the processed map at iteration i and the final it- 



eration, estimated as 2y C^Cb using auto-spectra. Here 

CI is the spectrum of the difference map between itera- 
tion i and the final iteration, number 1000. We define the 
convergence ratio Tc as this change in power given as a 
fraction of the uncertainty in the power, (T(Cb), and find 
it to be sufficiently small (less than 0.5) by iteration 500 
at all scales. The cross-correlation with WMAP suggests 
that the maps are well converged down to ^ ^ 200. We 
do not divide the angular power spectrum by a transfer 
function at any value of £. 

We test the isotropy of the power spectrum by estimat- 
ing the power as a function of phase 9 = ta.n~^ {£y fix). 
We compute the inverse-noise-weighted two-dimensional 
pseudo spectrum coadded over map subsets and patches. 
The mean cross-power pseudo spectrum is shown in Fig- 
ure [HI indicating the region masked at \£x\ < 270. The 
spectrum is symmetric for £ to as it is for any real- 
valued maps. To quantify any anisotropy, the power av- 
eraged over all multipoles in the range 500 < £ < 8000 
is computed in wedges of AO = 20°. It is found to be 
consistent with an isotropic 2D spectrum. 



TABLE 1 

Anisotropy power in bands of width Al = 300. 
Temperatures are in CMB units. 



Central 


£i£+l)Cb/2TT 


ACDM 


SZ 


Source 


Total 


4 




model'^ 


model'* 


model° 




750 


2317 ±98 


2106.9 


3.6 


0.8 


2111.2 


1050 


1313 ±42 


1204.2 


4.6 


1.5 


1210.3 


1350 


882 ± 20 


777.3 


5.3 


2.4 


785.0 


1650 


450 ± 22 


431.9 


5.8 


3.6 


441.3 


1950 


284 ± 14 


248.3 


6.1 


5.0 


259.5 


2250 


153 ± 12 


133.8 


6.3 


6.7 


146.8 


2550 


77 ± 9 


72.3 


6.5 


8.6 


87.4 


2850 


65 ± 9 


38.0 


6.6 


10.8 


55.3 


3150 


43 ±8 


19.8 


6.7 


13.1 


39.6 


3450 


29 ± 7 


10.3 


6.7 


15.8 


32.7 


3750 


50 ± 7 


5.3 


6.7 


18.6 


30.6 


4050 


47 ±12 


2.9 


6.6 


21.7 


31.2 


4350 


35 ±14 


1.6 


6.6 


25.0 


33.2 


4650 


34 ±19 


0.9 


6.5 


28.6 


36.0 


4950 


27 ±22 


0.6 


6.4 


32.4 


39.4 


5250 


77 ±33 


0.4 


6.3 


36.5 


43.1 


5550 


38 ±25 


0.3 


6.2 


40.7 


47.3 


5850 


42 ± 30 


0.2 


6.1 


45.3 


51.6 


6150 


26 ± 33 


0.2 


6.0 


50.0 


56.2 


6450 


103 ± 41 


0.1 


5.9 


55.0 


61.1 


6750 


-30 ± 68 


0.1 


5.8 


60.3 


66.2 


7050 


209 ± 115 


0.1 


5.7 


65.7 


71.6 


7350 


60 ± 139 


0.1 


5.6 


71.4 


77.1 


7650 


-38 ±213 


0.1 


5.5 


77.4 


83.0 


7950 


-63 ± 349 


0.1 


5.4 


83.6 


89.1 



''There is negligible covariance between bins. The maximum- 
likelihood fit agrees best with ACDM if the Cj, data given here 
arc multiplied by 0.91. 

For comparison with the SPT results of ILueker eFaD l(2009l 1. 
we also compute the spectrum with wider bins. For five bins of 
width Ai = 400 from 3000 < £ < 5000, we find 39 ± 6, 41 ± 7, 
49 ± 10, 36 ± 12, and 35 ± 13. For four bins of width Al = 900 
from 5000 < £ < 8600, we find 56 ± 19, 45 ± 25, 140 ± 70, and 
80 ± 250 loK^ . The ACT data are also shown in Figurc[T] with bins 
of width Al = 600 starting at £ = 4200 - 4800; the bandpowers are 
given by 31±11, 39±18, 36±19, 61±26, 29±58, and 28±116 [oK^. 

'^This and all "model" columns are in the same units as the data: 
thermodynamic [oK^. 

Assumes the best-fit value of ylgz = 0.63. 

•^Assumes the best-fit value of = 11.9 and uncorrclated Poisson 
sources. 
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Fig. 4. — The observed power spectrum in bandpowers at 150GHz from ACT observations (points with error bars). At large angular 
scales there is good agreement with the lensed ACDM model of the primary CMB (light blue curve shown for I < 4700). The of the 
model is 7.1 for 4 ACT data points in the range 600 < £ < 1800. The best-fitting model to the full dataset is shown (dark blue, the highest 
curve for all £ > 2500). The complete model includes the primary CMB model plus both a Poisson power from point sources and SZ power 
from clusters; both additional components have been allowed to vary. The complete model has been smoothed by convolution with a boxcar 
window function of width A£ = 300; the primary CMB model has not been smoothed. The narrower, gold band shows the marginalized 
95% CL limits on the Poisson amplitude, while the curve indicates the best-fit amplitude Ap = 11.9 |aK^. The wider pink band shows the 
95% CL upper bound on the SZ amplitude, Agz < 1.63; the dark curve inside it shows the best-fit value of Agz = 0.63. The Poisson and 
SZ power arc consistent with higher frequency observations and with ACDM predictions. The fitting procedure is described in Section [7] 



We test that the signals in separate data subsets are 
consistent by taking the cross-spectrum of two difference 
maps formed from the temperature maps, T*, of the four 
data subsets via: 

T^^{n) = [T^{h) -T^{h)]/2 (11) 
T34(n) = [T3(n) -r'*(n)]/2. (12) 

(The data subsets are described in Section 15.10 The 
difference maps are expected to contain noise but no 
residual signal. We estimate the cross-spectrum of the 

difference maps, Cb = (t'^'^T^^^ using the methods de- 
scribed in Section [5] for i\/ = 2 data segments. The two 
other permutations of the data, Ci, — ^T-'^^T^^^, and 

^14^23^ are also tested. The three difference spectra, 
shown in Figure [3 are consistent with no signal. 

7. CONSTRAINTS ON SZ AND IR EMISSION 

We perform a simple analysis of the power spectrum 
to quantify the combined contribution from dusty galax- 
ies and radio sources, and the level of SZ emission. 
We assume a ACDM cosmology with Icnsing of the 
CMB and parameters from the five-year WMAP anal- 



ysis combined with B AO and supernovae measurements 
(jKomatsu et al.l [20091 ). We defer a full investigation of 
cosmological parameter constraints until we improve the 
absolute calibration and better account for astrophysical 
foregrounds. 
Our model for the power at 148 GHz is 

where Bf^^ is the lensed primary CMB power spectrum; 
Bf^ is a template spectrum corresponding to a prediction 
for the SZ emission in a model with erg = 0.8 at 148 GHz; 
Ap quantifies the Poisson point source power, required 
to be positive; and Bf" corresponds to correlated point 
source power from clustered galaxies. The SZ template 
we use includes the correlated thermal and kinetic SZ 
(kSZ) effect derived from numerical sirn ulations and is 
described in detail in lSehgal et al.l (j2010t) . Its amplitude 
is assumed to scale with erg as the seventh power, such 
that Asz = ((Tg^/0.8)^ for fixed baryon density. This ac- 
counts approximately for the frequency-dependent com- 
bination of the thermal SZ component scaling as the 7.5 
to 8th power of erg and the sub-dominant kSZ scaling as 
the fifth power. The expected point source power Cp is 
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Fig. 5. — Convergence of the maps as a funetion of iteration. 
The mapmaking algorithm converges well by iteration 500. For 
signal-only simulations {top), the amplitude of fluctuations in the 
difference between the processed output map at iteration i and the 
input map (denoted (C^)^^^), is less than 1% of the amplitude of 
fluctuations in the input map {\/Cb) at all scales by i = 500. Itera- 
tions i = 5, 20, 100, and 500 are shown. For simulations with noise 
and for the data {middle and bottom, respectively), convergence is 
tested by estimating the maximum change in power between the 
processed map at iteration i and iteration 1000, as a fraction of the 
uncertainty in the power in the final map. For iteration 500 this 
fraction rc (described in Section ^ is sufficiently small, less than 
0.5 at all scales. 





I 



BOOO -4000 -2000 



2000 4000 6000 



Fig. 6. — The estimated two-dimensional power spectrum C^, 
multiplied by a factor of £ to emphasize the angular, rather than 
the radial, variation. The power is consistent with being isotropic, 
when divided into wedges of A9 = 20°. The vertical lines indi- 
cate the narrow region \£x\ < 270, where excess power from scan- 
synchronous signals contaminates the power spectrum. This re- 
gion is not used for the power spectrum analysis. The regions near 
\£y\) ~ (4000, 4000) are more noisy but not biased. A power 
spectrum computed without these regions is consistent with the 
one we present. 



given by 



Cr) 



ab 



(14) 



an integral over all sources with flux S up to some maxi- 
mum flux S'cut , where g{v) = {c^/2kv'^)x [(e^ - l)^/x^e^] 



^^g'l'» »»» i» 
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+ 
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o ^ i»» » ' » » »»»»^' » »a ^ $$j| 
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Multipole / 

Fig. 7. — Cross-spectra of difference maps formed from four 
data subsets. In the top panel, two maps T^^ = — and 
ji34 ^ j'S _rp4 j^j-g formed; both are expected not to contain signal. 
The cross-spectrum is consistent with no signal. The lower two 
panels show the same cross-spectra for the other two permutations 
of the four data subsets and are also consistent with no signal (x^ is 
25.2, 27.1, and 28.4 in the three panels with 25 degrees of freedom). 

converts flux density to thermodynamic temperature, 
and X = hv/kTcMB- Then Ap is the binned £(£-|-l)Cp/27r 
Poisson power at pivot £o = 3000. In thermodynamic 
I^K^, Cp = 0.698 X IQ-'^Ap. The conversion to Jy^sr"! 
at 148 GHz is CpiJy^sr-i] = l.hhCp[lQ-^\xK^ si]. 

The infrared sources are expected to be clustered. 
At small angles galaxies clu ster with typi cal correlation 
function C{e) oc O'^-^ (e.g.. IPeeblesl[T980l) . which would 
give Gi oc on non-linear scales. Motivated by this, 
wc first adopt a simple template for the correlated power, 



B'i 



I 



3000 



(15) 



and fit for the amplitude A^. Note that on larger 
scales, I < 300, B'i" " is expected to flatt en and gradu- 
ally turn over (e.s:.. IScott fc Whit^ I1999I) . but at these 
scales at 148 GHz the CMB dominates. Models sug- 
gest that the power from source clustering is less than 
the Poisson component at scales smaller than i = 2000 
(IScott fc White! fl99l INeerello et all [2007t IRiehi et alJ 
I2008D . This is consistent with observations at 600 GHz by 
BLAST (|Dcvhn et al. 2009; Vier o et al... 20091) . We there- 
fore impose a prior that the correlated power be less than 
the Poisson power at ^ = 3000, i.e., < Ac < Ap. This 
model is likely too simplistic, and t he correlated power 
may have an alternative shape fe.g.. lSehga.1 et al.l 120101 
or the halo model considered in I Viero et al.ll2009( ). If the 
power from dusty galaxies is instead better described by 
a linear matter power spectrum scaled by a bias factor, 
Bg°" becomes almost degenerate with the SZ component 
Sf^ (see, e.g., iHaU et al1l2009D . Since it is not possi- 
ble to separate these components with data at a single 
frequency, our limit on the SZ contribution at 148 GHz 
should be considered as the upper limit on the sum of 
the SZ power and a degenerate correlated source power. 
Multi-frequency data will enable us to investigate the 
shape and amplitude of the correlated component more 
fully 

Given this model, the likelihood of the data is given by 



2ln^ = (6^ 



06)+ In dots, (16) 



with covariance matrix S defined in Equation |8l The 
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Fig. 8. — Probability distributions for the amplitude Ap of an IR Poisson point source contribution of the form Bi = £{£ + 1)Ch/2it = 
Ap{£/3000)^, and of Agz, the amplitude of the SZ effect relative to that of a model with erg = 0.80. Each panel shows two models for the 
point source power spectrum. In one case, we assume the point sources are uncorrelatod Poisson-distributed sources. In the second case, 
wo marginalize over the amplitude of a correlated term scaling as Bt oc £. Left: The two-dimensional distribution of Ap and Agz- The 
filled (blue) regions assume uncorrelatod sources; the unfilled (rod) regions allow for correlated sources. Center: The distribution of Ap 
marginalized over Asz- In the center and right panels, the curve peaking at higher amplitude (blue) assumes uncorrelatod sources. The 
point source power is consistent with SCUBA and BLAST data at higher frequencies. Right: The distribution of Asz marginalized over 
Ap. The shaded region shows the 95% CL limits on erg inferred from WMAP data combined with distance priors. 



theoretical spectrum is computed from the model 
using Bf^ = WbiBf^, where wm is an approximate form 
of the bandpower window function in band h. We 
marginalize over the calibration uncer tainty analytically 
(IGanea et alJfToM IBridle et alJ[200l . 

The uncertainty on the shape of the window function is 
small, of order 1.5% for the 148 GHz band, compared to 
12% overall calibration uncertainty in power. The win- 
dow uncertainty is therefore neglected throughout this 
analysis. We verify that this approximation is valid by in- 
cluding the window function uncertainty in the likelihood 
calculation. The beam Legendre transform is first ex- 
panded in orthogonal basis functions, and the uncertain- 
ties on the basis function coefficients are used to derive 
the window function covariance matrix. The covariance 
is dominated by a small number of modes, so a singu- 
lar value decomposition is taken. The ten largest modes 
are included in the likeli hood, following the m ethod de- 
scribed in Appendix A of lHinshaw et al.l ()2007[) . We find 
that including the window function uncertainty has only 
a negligible effect on parameter estimates (< O.OSct). The 
ACT/MB AC beam measurements and th e orthogonal 
functi on expansion are described in detail in lHincks et al.l 
(|2009l) . 

We confirm that the expected values of parameters 
^sz, and Ac are recovered from maps of the ACT 
simulated sky. We tested simulations with and with- 
out noise in the map, and simulations with realistic 
timestream noise run through the mapper. 

In the range 600 < £ < 1800, where foreground emis- 
sion and secondary effects are sub-dominant, the data are 
consistent with the lensed ACDM model alone, with 
= 7.1 for four degrees of freedom (dof). If we rescale the 
maps to check consistency with the model (multiplying 
temperatures by 0.96, a 0.7cr change in the calibration), 
then the value becomes 3.0 with 3 dof. We do not 
rescale the maps in the analysis that follows, though the 
scale factor is allowed to vary when we marginalized over 
uncertainties. 



Using the full range 600 < £ < 8100, we find marginal- 
ized constraints on Ap ~ 11.2 ± 3.3 [a.K^ (thus, Cp = 
(0.78 ±0.23) xIO-VK^) and Asz < 1-63 (95% CL). 
The minimum = 27.0 for 23 dof. Assuming the scal- 
ing of Asz as {a§^y , this implies an upper limit of 
crf^ < 0.86 (95% CL). The one- and two-dimensional 
distributions are shown in Figure |S1 with limits given 
in Table [21 Marginalizing over the possible SZ power, 
ACT detects a residual point source component at 3(t 
{Sx^ = 10). Marginalizing over a correlated term with 
Ac < Ap gives Ap = 9.7 ± 2.8 [xK^, and af^ < 0.84 with 
= 26.7 for 22 dof. Because Ac is forced to be pos- 
itive, its inclusion has the effect of lowering the limits 
on Ap and Asz- Nevertheless, we do not find evidence 
for a correlated component with current sensitivity lev- 
els. The estimated parameters vary by less than 0.6a 
when the minimum angular scale is varied in the range 
5000 < £max < 8000, or the SZ template is replaced by 
the spectrum of lKomatsu fc SeliakI (1200l . which IS ap- 
proximately 15% lower than our template in the relevant 
range of £ from 1000-5000. We note t hat there are mul- 
tiple models for pred icting Asz (e.z.. IBond et al.l 120051 : 
iKravtsov et al.|[2005l) and that the relation between Asz 
and (Ts is an active area of research. 

We also co mbine the ACT spec trum with the WMAP 
5-year data (jPunklev et al.l [2009D to constrain the six- 
parameter ACDM model (defined by the baryon den- 
sity, cold dark matter density, cosmological constant, 
optical depth to reionization, and the amplitude and 
scale dependence of primordial fiuctuations at fc = 0.002 
Mpc^^). We model the SZ and point source contribution 
using Equation 1131 neglecting a correlated component. 
We find similar results for the point source and SZ am- 
plitude in this extended model, with Ap = 11.5±3.2, and 
Asz < 1.66 (95% CL). The ACDM parameters are con- 
sistent with WMAP alone, with lOOfJ^/i^ = 2.27 ± 0.06, 
flch^ = 0.111 ± 0.006, = 0.738 ± 0.030, = 
0.964 ± 0.014, T = 0.086 ± 0.017, and 10^ A, = 2.4 ± 0.1. 
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TABLE 2 

Constraints on secondary anisotropies and extra-galactic foreground emission at 150 GHz. 





Asz 




Ap{l = 3000) 


Cp 


Cp 












10- 


Jy2 sr-1 




Poisson point sources 


< 1.63 


< 0.86 


11.2 ± 3.3 


0.78 ±0.23 


1.22 ± 0.36 


27.0/23 


Poisson + correlated point sources 


< 1.36 


< 0.84 


9.7± 2.8 


0.68 ±0.20 


1.05 ± 0.30 


26.7/22 



7.1. Comparison to other point source observations 

The residual source level, Cp ^ (0.78 ± 0.23) 
xlO-^li.K^, combines power from radio sources and 
dusty galaxies that were not removed by the mask. 
For the three radio source models discussed in Sec- 
tion m we expect re sidual power of Ap = 6.4 \xK^ 
(|Toffolatti et al.l 119981, a f ter re scaling by 0.4), Ap = 
4.1 uK^ (IDe Zotti et al.l [2005h . and Ap = 7 \iK^ 
(ISehgal et al.ll201fln . These models' amplitudes corre- 
spond to C^""^'" = 0.43, 0.29, and 0.49xlO-VK^ re- 
spectively. The correction for those few sources below 
20mJy that we mask can be neglected; it is smaller 
than the spread among the models. Given the uncer- 
tainties in the models, wc subtract the typical model 
from the total source level and infer that the component 
from residual IR sources lies in the approximate range 
0.2 < Cj]^ X 10^ < 1 M-K^. 

The ACT result on total point source power is similar 
to those of APEX at 150 GHz, whic h finds Cp = l.lto j 
xIO^VK^ (iReichardt etldl l2009al) , of ACB AR {Cp = 
2.7t2;6 x10-5|^K2, IReichardt etHI I2009b^, an d SPT 
{Cp = (0.74 ± 0.06) X 10-5 uK'^. IHaU et all 120091) . This 
last measurement employs a lower flux level for remov- 
ing discrete sources: 6.4 mJy versus our cut at 20mJy. 
While this difference means that the Cp presented will 
contain roughly three times as much power from radio 
sources as the SPT measurement, we nevertheless find 
that the total residual power due to point sources is con- 
sistent between the two results, given the conservative 
assumption that at least one-quarter of the point source 
power observed by ACT is due to dusty galaxies rather 
than to radio sources. 

The IR source models in the literature for 148 GHz 
make a range of pre dictions. For e x ampl e, for sources 
less than 20mJy the iLagache et ahl ( 20041 ) model gives 



Cp = 40 xlO-VK^ the EegreHoeraD (200?) model 
gives Cp = 1.2 xlO-5 ^K^, and the[Schgal ct al. (20lJ) 
model gives Cp = 17.5 xlO-^ i xK^. T he best agreement 
comes with the iNegrello et al.l ()2007D model. 

The emission from dusty galaxies in the rest frame can 
be modeled as So{i^) fx u^By{T), with emissivity index 
/3 ~ 1.5 and Planck function B^{T). The effective in- 
dex a, where S{v) oc v'^ {Ce oc i/^"), accounts for the 
redshift of the sources, the intrinsic temperature varia- 
tion T, and the index /3. SCUBA has observed emis- 
sion at 850 |i.m (353 GHz), where we expe ct a similar 
popu l ation of galaxies to contribute (e.g., iGreve et al.l 
12004 120081) . Using SCUBA galaxy numbe r counts and 
a model for dN/ dS, iScott fc Whitl ()1999[ ) estimate the 
Poisson power to be Cp = 730x 10- VK^(190 Jy^sr-^) 
for S'cut ~ 50 mJy. Combined with residual IR source 



level observed by ACT, this implies an effective spectral 
index of ai5o-35o between 2.6-3.3. This is consistent 
with emission from dusty star-burst galaxies at high red- 
shift a nd in line with pred i ctions bv I White &: Majumdarl 
(|2004| ) and INegrello et al.l (|2007| ). It also agrees with the 
a = 2.6 ± 0.6 index inferred from source fluxes mea- 
sured with MAMB O (1.2 mm) and SCUBA (850 ^m) 
(jCreve et al.l I2004D , and with a = 2.3 measured from 
AZTEC (1.1mm) and SCUBA (jChapin et al.|[2009l ). Ob- 
servations by BLAST at 500 \xm (600 GHz) have Poisson 
power (2.7±0.2) x 10^ Jy^sr-^ (IViero et al.ll2009l) . Com- 
bining this with the ACT data leads to an estimate of 
the effective index 2.7 < aiso-eoo ^ 3.6, consistent with 
findings by APEX. The consistency with ai5o-35o sug- 
gests that similar populations are being probed at these 
frequencies, although BLAST is sensitive to a lower red- 
shift range than ACT. 

7.2. Comparison to other SZ observations 

The ACT constraints on the SZ power indicate an am- 
plitude of fluctuations trf^ < 0.86 (95% CL). This re- 
sult is consistent with estimates that combine the pri- 
mordial CM B anisotropy with di stance measures, cts = 
0.81 ± 0.03 (|Komatsu et al.ll2009| ). and improves on SZ- 
inferred limits at 150 GHz from Bolocam (cr|^ < 1.57) 
(jSavers et all 120091), Boom erang (crp < 1.14 at 95% 
CD (iVeneziani et al.l 120091) . and APEX {af- < 1.18) 
(jReichardt et al.ll2009a|) . We do not see evidence for an 
excess of SZ power, in contrast to lower frequency ob- 
servations by CBI at 30 GHz which prefer a value 2.5a 
higher than the co ncor dance value (crf^ = 0.922 ± 0.047, 
ISievers et an[2009l) . 

The ACT result is also consistent with the recently 
reported erf ^ = 0.773 ± 0. 025 (Arz = 0.42 ± 0.21) from 
the South Pole Telescope (jLueker et al.ll2009f ). The SPT 
team prefers a form for the correlated point sources that 
is covariant with the SZ template in the £ — 3000 range 
(jHall et al.l l2009f ). If such a form is correct, then the 
Asz we report should be interpreted as an upper limit 
on correlated point sources plus the SZ effect. Analysis 
of the ACT'S 218 and 277 GHz data wih shed light on 
possible forms of the correlated component. 

The ACT results on erg are also consistent with several 
recent studies based on the analysis of ROSAT X-ray 
flux-selected clusters. We note that these X-ray clus- 
ter studies themselves are consistent with measures of 
(Tg from richness or weak-lensing selected cluster sam- 
ples (see references in the three articles cited here). 
IHenrv et al.l ([20091) find crg(17™/0.32)0-30 = 0.86 ± 0.04 
(for ftm < 0.32) using clus ter gas temperatures rn easured 
with the ASCA sateUite. IVikhlinin eFall (|2009( ) obtain 
crg(rJ™/0.25)° ''^ = 0.813±0.013 (stat)±0.024 (sys) using 
temperatures derived from Chandra observations. Simi- 
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larly, iMantz et alj (l2009f ) find that in spatially flat mod- 
els with a constant dark energy equation of state, ROSAT 
X-ray flux-selected clusters yield ilm = 0.23 ± 0.04, 
0-8 =0.82 ±0.05. 

After this article was completed, WMAP 7-year mea- 
surement s of the brightness of Mars and Uranus were 
released (jWeiland et al.ll2010( ). They suggest that Tjj « 
107 K, which is ^ 4% dimmer than the value used here. 
If the new Uranus temperature is adopted, then all ab- 
solute temperatures and source brightnesses in this work 
would be reduced by 4%, and all Ce values by 8%. 

8. CONCLUSIONS 

We have presented the first map and the power spec- 
trum of the CMB sky made using data from the Atacama 
Cosmology Telescope at 148 GHz. With this map we can 
compare to WMAP at degree angular scales and mea- 
sure point sources with a resolution of 0?02. With an 
unbiased estimator, we extract the power spectrum, C'i, 
over a range of power exceeding 10^. 

We have interpreted the spectrum with a simple model 
composed of the primary CMB, a possible SZ contribu- 
tion, and uncorrelated point sources. This analysis pro- 
vides a new upper bound on the SZ signal from clusters 
(crf^ < 0.86 at 95% CL, though this is subject to un- 
certainty in the SZ models). A coordinated program of 
X-ray, optical, infrared, and millimeter-wavelength ob- 
servations of the largest SZ clusters is underway. 

These high angular resolution measurements probe 
the microwave power spectrum out to arcminute scales. 
Above i ~ 2500, the spectrum is sensitive to non-linear 
processes such as the formation of galaxy clusters and 
dusty galaxies. On the low-^ end, the spectrum mea- 
sures the Silk damping tail of the CMB which can be 
computed using linear perturbation theory as applied to 
the primordial plasma. It is clear that to understand the 
i > 1000 end of the primary CMB, and thus to improve 
significantly on measurements of the scalar spectral index 
and its running, source modeling will be required. Future 
analyses of the ACT data will include the two higher- 
frequency channels and additional sky coverage. In an- 
other approach, one can measure the high-^ E-modes, 
because the polarized CMB to foreground ratio is ex- 
pected to be higher than that for the temperature. We 
are pursuing both programs, as a polarization-sensitive 
camera is currently under development. 
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